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Abstract. It is shown that symmetric configurations of fuzzy spin direction detectors 
£\j ■ generate, through quantum jumps, IFS fractals on the sphere S* 2 . The IFS fractals can 

be also interpreted as resulting from applications of Lorentz boosts to the projective 
light cone. 
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1. Introduction 



Recently, Numerical Algorithms for quantum jumps were discovered. First they were 
introduced in quantum optics as convenient and effective tools for numerical simulations 
of the Liouville equation [1, 2, 3, 4]. A short history and more references can be found 
in Ref. [5]. Ph. Blanchard and A. Jadczyk, in a series of papers on EEQT - Event 
^ Enhanced Quantum Theory (cf. [6, 7, 8, 9, 10, 11] and references therein) - developed a 

new approach that unifies what John von Neumann called U- and R-processes [12] into 
a single piecewise deterministic process (PDP) where continuous evolution of a system 
is cyclically interrupted by discontinuous "jumps". 

Normally one would expect that time evolution of a physical system is described by 
a differential equation. That is how laws of physics are usually expressed. Here however 
we have a surprise: in EEQT a history of an individual quantum system, coupled to a 
monitoring device, as in every real world experiment, is described by a process or an 
algorithm rather than by a differential equation. The PDP is similar to those studied 
in the science of economics, where periods of smooth fluctuations are interrupted by 
market crashes [13, 14]. A quantum jump is what corresponds to a market crash - a 
discontinuity, a 'catastrophe'. But discontinuities and catastrophes have their own laws, 
and here comes the concept of a piecewise deterministic Markov process^ - a PDP. 

PDPs are the simple and elegant ways to describe the world in terms of cyclic, 
rather than linear, time; that is the world of cyclically, though somewhat irregularly, 

f The Markovian property is not really important. 
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recurring catastrophes. In the cases studied in the present paper the catastrophes come 
from the coupling of a quantum system to a system of two-state "detectors". In this 
case the catastrophes are not really catastrophes for the detectors - detectors just flip, 
which is exactly what detectors do for living. But these flips bring catastrophes for the 
quantum system, because with each flip of the detector, with each "event", as we call it, 
the quantum system state vector breaks its continuous evolution, and instantaneously 
jumps to a different state - it "rejuvenates" and it starts another cycle of a peaceful, 
continuous evolution - till the next catastrophe. 

Note. The term instantaneously in the last sentence may suggest that our formalism 
is incompatible with Einstein's relativity. That this is not so has been demonstrated 
in Ref. [10] (for a somewhat different approach cf also [15, 16, 17]). The point is that 
for a relativistic theory the role of time is being played by the Fock-Schwinger "proper 
time" as a Floquet variable. The reader should bear in mind that the EEQT algorithm 
is explicitly nonlocal: to simulate a history of an individual system integrations over 
entire space (or space-time) are needed. 

The EEQT algorithm generating quantum jumps is similar in its nature to a 
nonlinear iterated function system (IFS) [18] (see also [19] and references therein) and, 
as such, it generically produces a chaotic dynamics for the coupled system. Here the 
probabilities assigned to the maps are derived from quantum transition probabilities 
and thus depend on the actual point, but such generalizations of the IFS's have been 
also studied (cf. [20] and references therein). In the present paper we describe the 
algorithm generating quantum fractals, that is self-similar patterns on the projective 
plane Pi(C) ~ S 2 , when a continuous in time "measurement" of several spin directions 
at once takes place. 

Note: The term quantum fractals has been used before by Casati et al. [21, 22] in a different 
context. 

Note: The operators for different spin directions do not commute, but this does not contradicts 
Heisenberg's uncertainty relations as these deal with statistical ensembles averages, while here we are 
describing an individual quantum system. In fact, realizing the chaotic behavior of a quantum state 
vector, when several noncommuting observables are being simultaneously monitored, can help us to 
understand the mechanisms of statistical uncertainties. 

As stated above the algorithm of EEQT describes a piecewise deterministic random 
process - periods of a smooth evolution interspersed with catastrophic "jumps." Of 
course, once we have individual description, we can also get the laws for statistical 
ensembles. Here we get a nice, linear, Liouville evolution equation for measures - as it 
is usual in studying chaotic dynamics. The fact that there is a unique PDP generating 
the Liouville equation in the framework of EEQT has been proven in Ref. [23]. 
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2. EEQT - Quantum Fractals 

2.1. Geometry 

We have published, as an OpenSource project [24], the algorithm implemented in Java 
that generates the Five Platonic Fractals - that is fractals generated by five most 
symmetric detector configurations. The algorithm generates self-similar patterns on 
a sphere of a unit radius. The points on the sphere represent (pure) states of the 
simplest quantum system - the spin 1/2 rotator. This spin 1/2 quantum system is 
coupled, continuously in time, to a finite number of symmetrically distributed spin- 
direction detectors. Thus the symmetry of the pattern reflects the symmetry of the 
detector directions distribution. Each spin direction is characterized by a vector n of 
unit length. Here we study the most symmetrical configurations, therefore we chose 
direction vectors nj pointing from the origin to the vertices of one of the five platonic 
solids. We consider the following five detectors configurations: 

(i) tetrahedron: 4 detectors along the directions n[i],i = 1, . . . , 4 

{{0, 0, 1.}, {a[17], 0, -a[3]}, {-a[6], a[12], -a[3]}, {-a[6], -a[12], -a[3]}} 

(ii) octahedron: 6 detectors along the directions n[i], % — 1, . . . , 6 
{{0, 0, 1.}, {1., 0, 0}, {0, 1., 0}, 

{-l.,0, 0}, {0, -1., 0}, {0, 0,-1.}} 

(iii) cube: 8 detectors along the directions n[i], % = 1, . . . , 8 

{{0, 0, 1.}, {a[17], 0, a[3]}, {-a[6], a[12], a[3]}, {-a[6], -a[12], a[3]}, 
{a[6], a[12], -a[3]}, {a[6], -a[12], -a[3]}, fa[17], 0, -a[3]}, {0, 0, -1.}} 

(iv) icosahedron: 12 detectors along the directions n[i], % = 1, . . . , 12 
{{0, 0, 1.}, {0.a[15], 0, a[5]}, {a[2], a[13], a[5]}, {-a[10], a[7], a[5]}, 
{-a[10], -a[7], a[5]}, {a[2], -a[13], a[5]}, {a[10], a[7], -a[5]}, 
{a[10], -a[7], -a[5]}, fa[2], a[13], -a[5]}, {-a[15], 0, -a[5]}, 
{-a[2],-a[13],-a[5]}, {0, 0, -1.}} 

(v) dodecahedron: 20 detectors along the directions n[i],i = 1, . . . , 20 
{{0, 0, 1.}, {a[9], 0, a[ll]}, {-a[3], a[8], a[ll]}, {-a[3], -a[8], a[ll]}, 
{a[ll], a[8], a[3]}, {a[ll], -a[8], a[3]}, {-a[14], a[4], a[3]}, 

{a[l], a[16], a[3]}, {a[l], -a[16], a[3]}, {-a[14], -a[4], a[3]}, 
{a[14], a[4], -a[3]}, {a[14], -a[4], -a[3]}, {-a[ll], a[8], -a[3]}, 
{-a[l], a[16], -a[3]}, {-a[l], -a[16], -a[3]}, {-a[ll], -a[8], -a[3]}, 
{a[3], a[8], -a[ll]}, {a[3], -a[8], -a[ll]}, {-a[9], 0, -a[ll]}, 
{0, 0, -1.}} 

where the array of real numbers a[i],i = 1, . . . , 17 is given in the following table. 
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Note: There is no deep reason why we chose these configurations - simplicity and beauty are the main 
factors here. Notice that, to enable easy zooming onto the attractor, we have chosen the orientations 
in such a way that in each case the North Pole, with coordinates (0,0,1), of the sphere is occupied by 
one of the vertices. 



2.2. The algorithm 

Here we describe the algorithm. Comments on its meaning and on the derivation can 
be found in the endnotes in Section 4 

2.2.1. The Hilbert space The quantum system is represented in a two-dimensional 
complex Hilbert space, which we realize as C 2 - the set of all column vectors 



(1) 

where a\ and a 2 are complex numbers, and with scalar product (a, b) defined by 

(a, b) = ai&i + a 2 b 2 (2) 
where the bar c stands for the complex conjugation of the complex number c. 

2.2.2. Spin directions We choose the Pauli matrices cr x ,a y ,a z to represent spin 
directions along x, y, z axes respectively. 

fo, 1\ /0, -i\ fl, 



i, oy z y \i, o y 3 2 vo, -i 

Together with the identity matrix 

— 0; !) 

they span the whole 2x2 complex matrix algebra. In computations it is important 
to make use of the fact that Pauli matrices (after multiplication by "-i") represent the 
quaternion algebra, that is: 

Jl J2 ^2 T 

a±a 2 = —o- 2 o- 1 = ia 3 , a 2 a 3 = ~cr 3 a 2 = i<7i, 0-30-! = —0\0 3 = \o 2 
To each direction n in space there is associated spin matrix 

/ \ . . f n 3 , m - in 2 \ 

cr(n) = n x a x + n 2 a 2 + n 3 a 3 = (4) 

\n 1 + in 2 , -n 3 J 
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satisfying automatically c(n) 2 = /, and with eigenvalues +1,-1. Vectors {^j an< ^ 

1^ are eigenvectors of <r 3 to eigenvalues +1 and —1 respectively and thus correspond 

to "North" and "South" spin orientations respectively. Let -P(n) denote the projection 
operator that projects onto eigenstate of <r(n) to the eigenvalue +1. Then P(n) is given 
by the formula: 

P(n) = i(/ + (7(n)). (5) 

Indeed, P(n) is Hermitian and has eigenvalues |(1±1) = 1 or - thus it is the orthogonal 
projection, and it projects onto the eigenstate of cr(n) with spin direction n. 

2.2.3. Fuzzy projections For each n let P(n, e) be the fuzzy projection opearator 
defined by the formula: 

P(n,e)= 1 -(I + ea(n)) (6) 

where < e < 1 ), or better: 1 — e, is a parameter that measures the "fuzziness." The 
extreme cases are not the very interesting ones: for e = we get the identity operator 
- maximal fuzziness and no information whatsoever, while for e = 1 we get the sharp 
projection -P(n) = P(n, e = 1). 

We restrict the range of the parameter e to the interval [0, 1] because only in this 
range P(n, e) is a positive operator. It is easy to see that this is so. Indeed, a Hermitian 
matrix is positive when its eigenvalues are positive, and the eigenvalues of P(n, e) are 
(1 ± e)/2, thus — 1 < e < 1. On the other hand negative e for n is the same as positive 
e for — n, thus we restrict the range of e to [0, 1]. 

It is the operators P(n, e) that will act on quantum states to implement "quantum 
jumps" whenever detectors "flip." 

The overall coefficient in the definition (5), chosen to be \ here, is not important 
because in applications each of the operators P(n, e) is multiplied by a coupling constant, 
and, in our case, when we are not interested in timing of the jumps, the value of the 
coupling constant plays no role. 



2.2.4- Jumps are implemented by fuzzy projections Let us now discuss the algebraic 
operation that is associated with each quantum jump. Suppose before the jump the 
state of the quantum system is described by a projection operator P(r), r being a unit 
vector on the sphere. That is, suppose, before the detector flip, the spin "has" direction 
r. Now, suppose the detector P(n, e) flips, and the spin right after the flip has some 
other direction, r'. What is the relation between r and r'? It is easy to see that the 
action of the operator P(n, e) on a quantum state vector is given, in terms of operators, 
by the formula: 

A(e, n, r)P(r') = P(n, e)P(r)P(n, e), (7) 
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where A(e, n, r) is a positive number. It is a simple (though somewhat lengthy) matrix 
computation that leads to the following result: 

A(e,n,r) = l±^±^ (8) 

(l-e 2 )r + 2e(l + e(n-r))n 

1 + e 2 + 2e(n • r) K ' 

where (n • r) denotes the scalar product n • r = nir\ + n 2 r 2 + n 3 r 3 . 

2.2.5. Transition probabilities Given the actual state r of the quantum system, and 
the configuration of the detectors {n[?],i = 1,2,..., TV} we compute probabilities 
p[i], Y2iLiP[i] = 1 f° r the i — th detector to flip. Then we select randomly, with the 
calculated probability distribution, the flipping detector, and we implement the jump 
by changing r to r' according to the formula (7), with n = n[i]. The probabilities p[i] 
are computed using the theory of piecewise deterministic Markov processes applied to 
the case of quantum measurements - as developed within EEQT. | According to EEQT 
the probabilities p[i] are given by the formula: 

p[i] = const • Tr (P(r)P(n[i\, e) 2 P(r)) (10) 

where const is the normalizing constant. Using cyclic permutation under the trace, as 
well as the fact that P(r) 2 = P(r) we find, taking trace of both sides of the formula (7), 
that p[i] are proportional to A(e, n[i],r) given by (8), thus 

1 + e 2 + 2e(n[i] • r) 
N{1 + e 2 



p[«i = Jr. " J • (") 



Note that, owing to the fact that XlfcLi n [^] = we have ^2f = iP[i] = 1, as it should be. 



3. The Five Platonic Fractals 



As we noted above, different values of e give different fuzziness. To produce 
representative pictures, one for each solid, we adjusted e so that, as a rule, the more 
vertices, the higher value of e - thus higher resolution of details. While rendering the 
pictures, for obtaining grayscale value for a given pixel we were using either the formula 
log(data + 1) or, to get more details at the peak values, even log(log(data + 1) + 1). 

3.1. Other Poly hedra 

In principle our algorithm should create quantum fractals for each of the regular 
polyhedra. The only restriction on the array of vectors n[i] is that they are all of unit 
length, and their sum is a zero vector. We added, for comparison with the Platonic 

X It is of interest that Born's probabilistic interpretation of quantum mechanics as well as the standard 
formula for quantum mechanical transition probabilities, can be derived in this way and there is no 
need of adding it as a separate postulate. 
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solids configurations, two additional simple yet regular figures: double tetrahedron 
and icosidodecahedron. Notice that tetrahedron is self-dual, while dodecahedron and 
icosahedron are dual to each other. Double tetrahedron array is obtained by combining 
n[i] with —n[i] - that is with the inverted configuration. 

Icosidodecahedron has particularly simple and elegant expression for its 30 vertices: 
they are of the form: (±1,0,0) and its cyclic permutations, and |(±1,0, ^r) an d its 
cyclic permutations, where = 1+ 2 V ^ = 1.61803. . . is the golden ratio. All of its edges 
are of length 0. For its 30 vertices e = 0.85 was needed to resolve the atrractor's fine 
structure. 

4. Notes 

4-0.1. Complex projective plane serves as a canvas Pure states of the quantum spin 
(we are discussing spin 1/2 here) are described by unit vectors in our Hilbert space C 2 . 
But, as it is standard in quantum theory, proportional vectors describe the same state 
- the overall phase of the vector has no physical significance. Therefore, in geometrical 
terms, the set of all pure states is nothing but the projective complex space Pi(C) 
which happens to be the same as the sphere S 2 . That is why the fractal pattern, in 
our case, is being drawn on a spherical canvas. There is, however, another possible 
interpretation of the same algorithm. It is well know that there is an intrinsic relation 
between Minkowski space and the space of 2 x 2 complex Hermitean matrices. § In 
coordinates the map is given by p = i— > = p^cr^, so that det(^) = p 2 = p^p^. In 
particular our projection operators -P(r), r 2 = 1 correspond to null directions. In other 
words our canvas, the sphere S 2 , can be also thought of as the projective light cone - 
the space of light directions. Quantum jumps would then correspond to sudden changes 
of directions of light or light-like entities. Indeed the operators P(n, e) are positive, 
thus proportional to Lorentz boosts. The formula (9) for jumps implemented by these 
operators should be compared with the formula for a Lorentz boost, with velocity (3 in 
the direction n: 

y° = x° cosh a + (xn) sinh a, y = x — [(xn) cosh a + x° sinh a]n, 

where (3 = tanha is the velocity. It follows that the jumps described by the equation 
(9) can indeed be interpreted as Lorentz boosts witht velocity (3 = 2e/(l + e 2 ). 

The formula (7) can be easily generalized for P(r) being a generic Hermitean matrix 
(thus representing a four-vector p of space- or time-like character as well. However there 
is no such generalization for the formula (10), so that the probabilities would have to 
be assigned equal, and a physical interpretation, even tentative one, is missing in such 
a case (cf Subsection 4.11 below). 

§ In fact, by using Cay ley transform, this relation identifies the space of unitary matrices with the 
compactificd Minkowski space. 
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4-1. Pure states as projection operators 

For a spin 1/2 quantum system pure quantum states are uniquely described by projection 
operators P(n), where n is a unit length direction vector, starting at the origin, and 
ending on one of the points of the unit sphere. All pure states are of this form. Indeed, 
every orthogonal projection, except of the two trivial ones: and /, are of the form 
P(n) for some n. 

4-2. Fuzzy projections 

Detecting the particle spin is somewhat different than detecting its position components. 
To measure the position of a particle we can use a photographic plate or a bubble 
chamber. In such simple mathematical model of a detector is obtained by 

associating with each active center of the detector a fuzzy, bell-type function with its 
half-width corresponding to the active region of the center ( for instance about 1/xm for 
an AgBr grain of a photographic emulsion). What would correspond to such a "fuzzy" 
projection operator in the case of a spin measurement? We do not have much of a choice 
here. Due to symmetry reasons there is only one formula possible, namely one given by 
Eq. (6). 

4-3. Importance of fuzziness 

It is importance to note that in our generalization of the projection postulate, as the 
result of a jump, not all of the old state is forgotten. The new states depends, to some 
degree, on the old state. Here EEQT differs in an essential way from the naive von 
Neumann projection postulate of quantum theory. The parameter e becomes important. 
If e = 1 - the case where P(n, e) = -P(n) is a projection operator - the new state, after 
the jump, is always the same, it does not matter what was the state before the jump. 
There is no memory of the previous state, no "learning" is possible, no "lesson" is taken. 
This kind of a "projection postulate" was rightly criticized in physical literature as being 
contradictory to the real world events, contradicting, for instance, the experiments when 
we take photographs of elementary particles tracks. But when e is just close to the value 
1, but smaller than 1, the contradiction disappears. This has been demonstrated in the 
cloud chamber model [9], where particles leave tracks much like in real life, and that 
happens because the multiplication operator by a Gaussian function does not kill the 
information about the momentum content of the original wave function. Notice that 
our fuzzy projection operators P(n, e) have the properties similar to those of Gaussian 
functions, namely|| 

P(n,e) 2 = ^P(n >r ^). (12) 

|| This can also be interpreted as Lorentz formula for addition of relativistic velocities- cf. Section 4.0.1 
below. 
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4-4- Geometrical meaning of the parameter e 

It is instructive to have a visual picture of the map r i— > r' of the sphere S 2 implemented 
by the operator P(n, e). To this end let us assume the vector n is pointing North, i.e. 
n = (0,0, 1). Then the result r' of applying the operator P(n, e) to a point r on the 
sphere is on the same longitude as the original point r, but its latitude 9 changes - it 
moves towards the North Pole along its meridian, the new latitude being given by the 
formula: 

(l-6 2 )cOs(fl)+26(l + 6COs(fl)) \ 

l + e 2 + 2ecos(fl) ) 1 ' 



45 90 135 180 




45 90 135 180 



& 



arccos 



Figure 1. The amount of shift 9 — 0'(6), as a function of 9 for different values of e. 

Remark: Here 9 is not exactly the "geographical latitude". It is zero at the "North 
Pole" ((73 = +1), 90 degrees at the equator, and 180 degrees at the "South Pole" 

(ffs = -l)- 

Each map P(n, e) maps the sphere onto itself in an injective way. For n = (0, 0, 1) 
the map is easy to picture. All points of the sphere move towards the North Pole 
along their meridians, except of the two fixed points: North and South Pole. All of 
the Northern hemisphere, and a strip below the equator, shrinks, while the other part, 
near the South Pole, stretches. The amount of stretching can be found by plotting the 
function 9 — 9' (9) - it has a maximum at 9 corresponding to z — — e. Thus the parameter 
e gets a simple interpretation: it is the value of z coordinate for which shrinking of 
meridians is replaced by stretching - an equilibrium point. This point is always on the 
southern hemisphere. For e close to zero, where the map is close to the identity map, 
the equilibrium point is close to the equator. Then, as e approaches the value of 1.0, 
corresponding to the sharp projection operator, the equilibrium latitude gets closer and 
closer to the South Pole. In the limit of e = 1 all of the sphere shrinks to the North 
Pole, only the South Pole remains where it was. 

4-5. Quantum Fractals and IFS 

Our algorithm is, in fact, a version of a nonlinear iterated function system (IFS). Such 
algorithms are known to produce complex geometrical structures by repeated application 
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of several non-commuting affine maps. The best known example is the Sierpinski 
triangle generated by random application of 3 x 3 matrices A[i], % = 1, 2, 3 to the vector: 
vo = (^o, 2/o,l) where A[i] are given by A[i] = ((0.5, 0, axi), (0, 0.5, ayi), (0, 0, 1)) and 
axi = 1.0, ayi = 1.0, ax2 = 1.0, ay2 = 0.5, 0x3 = 0.5, ay 3 = 1.0. (Our 3x3 matrices 
encode affine transformations - usually separated into a 2 x 2 matrix and a translation 
vector.) At each step one of the three transformations A[i],i = 1,2,3 is selected with 
probability p[i] = 1/3. After each transformation the transformed vector is plotted on 
the (x, y) plane. Theoretical papers on IFSs usually assume that the system is hyperbolic 
that is that each transformation is a contraction, i.e. the distances between points get 
smaller and smaller. It was shown in [20] that this assumption can be essentially relaxed 
when transformations are non-linear and act on a compact space - as is in the case of 
quantum fractals we are dealing with. But it is not known whether the results of [20] 
apply in our case. 

4-6. Continuous part of the evolution 

In our discussion of quantum fractals, we neglect completely the continuous part of 
the time evolution. Here we are not interested in "when" jumps happen. We are only 
interested in the final pattern produced by a long sequence of jumps. Normally spin 
interacts with a magnetic field (if present), which causes our spin sphere to rotate around 
the direction of the magnetic field vector. Such a rotation, between jumps, would smear 
out our pattern. To study the pattern we neglect the rotation part. Once we discard 
the continuous evolution, the timing of jumps does not really matter, so we neglect this 
part of EEQT algorithm (timing is very important in simulations of particle detectors, 
arrival and tunneling times. Here, we simulate jumps as fast as our computer can crunch 
the numbers. 

4-7. Nonunitarity 

There is one important comment that applies here. Even if we neglect the continuous 
time evolution due to magnetic field, the very presence of the detectors causes a non- 
unitary time evolution of the pure state of the quantum system. This evolution is also 
called, by some physicists (Dicke, Elitzur, Vaidman) , "interaction-free". We do not 
want to enter into this subject here, except for one remark: for symmetric geometric 
configurations that we are considering here, the EEQT algorithm implies that this, 
continuous in time, non-unitary evolution can be neglected as well. In fact, it follows 
from the EEQT model that the "interaction-free" or, as we call it, "binamical part" 
of the evolution is determined by the generator —re A, where A = P(n[i], e) 2 . In our 
case, when our detectors are symmetrically placed, so that J^Li n W = 0, the formula 
(12) implies A = N ^ e - , thus the "binamical" part is just decreasing the norm of the 
state vector, while leaving its direction unchanged. Thus it does not affect the geometric 
pattern of jumps (it is responsible for the mean frequency of jumps, but here timing is 
not important). For a recent review of EEQT, cf. [28] 
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4-8. Quantum characteristic exponent 

Averaging our nonlinear PDP over individual histories one gets a linear Liouville 
equation for the density matrix of the total system. Tracing over the classical subsystem 
is, in our case, easily performed and we then get: 

N N 

p = «(£ P(n\i], e) P P(n\i],e) - -{£ P(n\i], e)\ p}) (14) 
i=\ i=i 
where {, } stands for anti-commutator, and k is a coupling constant. For p written as 
p = \(I + c(m)), m 2 = m\ + m| + m\ < 1 after some calculations we get a very simple 
time evolution: m(t) = exp (— iV/te 2 t/3) m(0). The quantum characteristic exponent, as 
defined in Ref [26], is thus ^-k - not a very useful quantity in our case. The Hausdorff 
dimension of the limit set, for the tetrahedral case, has been numerically estimated in 
Ref. [27] and shown to decrease from 1.44 to 0.49 while e increases from 0.75 to 0.95. We 
hope that by publishing the generating algorithm we will create interest in confirming 
these as well as obtaining new results in this field. 

4-9. How to measure the wave function itself 

The fractal patterns are produced by a jumping point on the space of pure quantum 
states - thus are not directly observable. That is why we consider our model as simply a 
toy to play with. But the model can be developed further on so as to predict observable 
effects. One way to do it is by adding another layer of detectors, densely spaced, that 
would detect the fractal pattern. Here we come to the famous question: how to measure 
the wave function itself. The idea as to how to model such a measurement within the 
EEQT formalism has been indicated in [8]. We hope to implement the appropriate 
algorithms at a later time. On the other hand, interpreting the patterns as resulting 
from Lorentz boosts applied to directions in real space, as discussed in Section 4.0.1, 
we may try to find similar fractal patterns in the distribution of Galaxies (compare for 
instance Figure (6) with recent paper on the topology of the universe [29]. 

4-10. Detectors 

A detector is represented by a two-state classical system. It can be in one of the two 
states, denoted and 1. The fact that it is "classical" means that its two states define 
two different superselection sectors that can be mixed statistically, but there are no 
observables that connect these sectors. We will assume that it can "flip" from to 1 or 
from 1 to when coupled to a quantum spin. Each flip represents an event; specifically: 
a detection event. The interpretation is that when the detector flips, the experimental 
question "is the spin oriented along the vector n?" gets an affirmative answer. Note that 
n and — n are two different experimental questions. They corresponds to two opposite 
spin directions. 

A realistic detector should also exhibit a relaxation time, that is, after each flip 
it should take some time before it is ready to flip again. We could easily model this 



Quantum Jumps, EEQT and the Five Platonic Fractals 



12 



phenomenon in our model, but here we are interested in patterns that are created, not 
in the timing of its appearance . 

4-11. Hyperbolicity 

In a recent paper [30] Lozinski, Slomczynski and Zyczkowski studied iterated function 
systems on the space of mixed states, when probabilities that are associated with maps 
are given independently of the maps. In section III of their paper they give a short 
discussion dealing with iterated function systems on the space of pure states as well. 
They start with the following definition of a (pure states) quantum iterated function 
system(QIFS): 

Definition (QIFS) Let ^ be a complex Hilbert space of dimension N. Let Vn be 
the space of one-dimensional subspaces of S) N . Given a unit vector <p £ f)Ni let P<j> be 
the orthogonal projection onto the subspace spanned by <fi.. Specify two sets of k linear 
invertible operators: 

• Vi : £jjv — > $)n (i — 1, • • • , k), which generates maps Fj : Vn — > Vn (i — 1, ■ ■ ■ ,k ) 
by h-» Vi0/||Vi0|| for any £ V N , and 

• Wi : S)n —> $3n (i — !,-••,&), forming an operational resolution of identity, 
Yli=i W-Wi = I, which generates probabilities pi : Vn — > [0, 1] (i = 1, . . . , k) 
byp,(0) :=||Wi(0)|| 2 . 

Comments 

(i) The authors define the system to be hyperbolic if the maps Fi are contractions 
with respect to the Fubini-Study distance d(P^,P^) = arccos {^/Tr{P ( f > P^)^ , i.e. 
there exists constants < Li < 1 such that d{Fi{P^) , Fi(P^)) < L i d(P(j) ) P^) 
for all 0, ip £ i^Ar. Then they state a proposition (Proposition 1 in [30]) that 
guarantees existence of an invariant measure for a hyperbolic system. It seems 
that the assumptions of this proposition can not be satisfied. It is well known [31] 
that a smooth injective map of a compact orientable manifold is automatically a 
bijection. On the other hand, iterating the map several times if necessary, the 
distance between any two image points can make smaller than any given positive 
number, and therefore, a fortiori less than the maximal distance between the points 
of Vn, which contradicts surjectivity.^f 

(ii) The authors of [30] assume that the probabilities pi are given independently of the 
mapping operators Vi. This is not the case in the EEQT scheme. In EEQT Vi need 
not be invertible, but the probabilities pi are determined by the V^-s automatically, 
and in such a way that whenever there is a danger of dividing by zero, the associated 
probability is automatically zero. 

If After informing the authors about this problem, they kindly replied that they saw it too, that their 
Proposition 1 has been deleted from the paper they submitted to Physical Review, and that they are 
going to replace the electronic version of the paper on the eprint server as well. 
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Figure 3. Quantum Octahedron, e = 0.58 




Figure 5. Quantum Icosahedron. e = 0.75 
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Figure 7. Quantum Double Tetrahedron, e = 0.7 




Figure 8. Quantum Icosidodecahedron. e = 0.85 



